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We report a measurement of the inclusive electron energy spectrum for semileptonic decays of 
B mesons in a data sample of 52 million T(45') BB decays collected with the BABAR detector 
at the PEP-II asymmetric-energy B-meson factory at SLAC. We determine the branching fraction, 
first, second, and third moments of the spectrum for lower cut-offs on the electron energy between 
0.6 and 1.5 GeV. We measure the partial branching fraction to be B{B — > Xeu,Ee > 0.6 GeV) = 
(10.36 ± 0.06(stat.) ± 0.23(sys.))%. 



PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 

The operator product expansion provides correc- 
tions to the relation between the semileptonic B de- 
cay rate and the magnitude of the Cabibbo-Kobayashi- 
Maskawa (CKM) [3| matrix element Vet in the free-quark 
model Q . The corrections are expressed in terms of non- 
perturbative quantities that can be extracted from mo- 
ments of inclusive distributions. We plan to use the preci- 
sion measurements of moments of the lepton energy spe_c- 
tra presented here and of hadron mass distributions 3] 
to determine those parameters and thereby to improve 
the determination of jV^^I 

In this paper, we present a new measurement of 
the inclusive electron energy spectrum from semilep- 
tonic B decays, averaged over charged and neutral B 
mesons produced at the T(4S') resonance. After cor- 
recting for charmless semileptonic decays, we derive 
from this spectrum several moments as a function of 
the minimum electron energy ranging from 0.6 GeV 
to 1.5 GeV, where lower endpoint is set by the lim- 
its of electron identification and prevalence of back- 
ground. In the B meson rest frame, we define Ri{Eo,fi) 
as J^^{Ee — iiy^{dr/dEe) dE^, and measure the first mo- 
ment Mi{Eo) ^ Ri{Eq,0)/Ro{Eo,0), the central mo- 
ments M„{Eo) = RniEo, Mi{Eo))/Ro{Eo,0) for n = 2, 3 
and the partial branching fraction B{Eq) = tb Ro{Eo,0), 
where tb is the average lifetime of charged and neutral 
B mesons. 

The measurements presented here are based on data 
collected by the BABAR detector at the PEP-II asym- 
metric e~^e~ storage ring; they correspond to an inte- 
grated luminosity of 47.4 fb"^ on the T(45') resonance 
and 9.1 fb~^ at an energy 40 MeV below the resonance 
(off- resonance) , measured in the electron-positron center 
of mass frame. Where background and efficiency cor- 
rections cannot be measured directly from data, we use 
a full simulation of the detector based on GEANT4 [^. 
In the following, all kinematic variables defined in the 
T(4S') rest frame will be annotated with an asterisk. 

This analysis is similar to the BABAR measurement 
of the semileptonic branching fraction 7] , including use 
of the same electron identification criteria, but super- 



sedes it by an order of magnitude in integrated lumi- 
nosity. We identify BB events by observing an elec- 
tron, ctag, with charge Q{etag) and a momentum of 
1.4 < p* < 2.3 GeV/c in the T(4S') rest frame. These 
electrons make up the tagged sample that is used as nor- 
malization for the branching fraction. Each electron egig 
with charge Q{esig) for which we require p* > 0.5GeV/c 
is assigned to the unlike-sign sample if the tagged sam- 
ple contains an electron with Q{etag) = —Qissig), and to 
the like-sign sample if Q{etag) = Q{&sig)- In events with- 
out B'^B'^ mixing, primary electrons from semileptonic B 
decays belong to the unlike-sign sample while secondary 
electrons contribute to the like-sign sample. 

Multi-hadron events are selected by either requiring a 
track multiplicity Nch > 5, or Nch = 4 plus at least two 
photon candidates with E^ > 80 MeV. Track pairs from 
converted photons are not included in Nch, but count as 
one photon. For further suppression of non-BB events 
we require the ratio of the Fox- Wolfram moments H2 / Hq 
to be less than 0.8. 

Electrons originating from the same B meson as the 
tagged electron typically have opposite charge and direc- 
tion. To reject them we require 

cosa* > 1.0 -p*(GeV/c) and cos a* > -0.2, (1) 

where a* is the angle between the two electrons. This re- 
quirement also excludes electron pairs from J/ijj — > e^e~ 
decays. To suppress background contributions from 
J/tp — e"'"e~ decays to the tagged sample, we require 
the invariant mass M^e of the tag electron, paired with 
any electron of opposite charge and cos a* < —0.2, to be 
outside the interval 2.9 < Mee < 3.15GeV/c^. Here the 
requirement on cos a* does not reduce the efficiency of 
this veto, but ensures that no signal electron satisfying 
Eq. ^is excluded from the unlike-sign sample. The effi- 
ciencies of these selection criteria are estimated by Monte 
Carlo (MC) simulation. 

Continuum background is subtracted from the tagged, 
like- and unlike-sign samples by scaling the off-resonance 
yields by the ratio of on- to off-resonance integrated lu- 
minosities, corrected for the energy dependence of the 
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continuum cross section. In the ofF-resonance sample, 
the momenta are scaled by the ratio of the on- and off- 
resonance energies. 

Electron spectra from photon conversions and Dalitz 
decays are extracted from data, taking into account the 
pair-reconstruction efficiencies from MC simulation. The 
relative uncertainty in these efficiencies is estimated to 
be 13% and 19% for conversion and Dalitz pairs, respec- 
tively. 

The misidentification rates for pions, kaons, and pro- 
tons are extracted from data control samples. They rise 
from 0.05% to 0.12% for pions and faU from 0.4% to 0.1% 
for kaons as p* increases from 0.5 to 2.5 GeV/c. The 
systematic errors are estimated from the control sample 
purities and from the uncertainties in the tt, K and p 
abundances. The resulting relative uncertainties are less 
than 40%. 

There is a small residual background in the sample of 
unlike-sign pairs originating from the same B meson and 
fulfilling the requirement on the opening angle a* from 
Eq. n It is estimated from a fit to the cos a* distribu- 
tion, separately for each 50- MeV/c-wide bin in p* . The 
distribution is flat for signal pairs, while for background 
pairs it is taken from MC simulation, with a maximum at 
cos a* = — 1 and gradually decreasing to at cos a* = 1. 
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FIG. 1: Measured momentum spectrum (points) and esti- 
mated backgrounds (histograms) for electron candidates in 
(a) the unlike-sign sample, and (b) the like-sign sample. 



Fig- n shows the electron momentum spectra and the 
background contributions discussed so far. Further back- 
grounds arise from decays of r leptons, charmed mesons 



produced in 6 — > ccs decays and J/^ or ipi^S) — > e~^e~ 
decays with only one detected e. We also need to correct 
for cases where the tagged electron does not originate 
from a semileptonic B decay. These backgrounds are ir- 
reducible, and their contributions to the three electron 
samples are estimated from MC simulations, using the 
ISGW2 model '^l to describe semileptonic D and Dg me- 
son decays. Assuming T{Ds Xev) — T{D Xev), 
we obtain BWs Xev) = (8.05 ± 0.66)%. Using 
0.84 ± 0.09 for the measured fraction of B ^ DsX 
decays where the Ds originates from fragmentation of 
the W Boson, and B{B D^X) = (10.5 ± 2.6)% [ij 
yields 6(5°'+ ^ D+ e+) = (0.71 ± 0.20)%. As- 
suming equal production rates of D and D* and us- 
ing B{B DD(*'>X) = (8.2 ± 1.3)% 113, we arrive at 
S(B°'+ ^ ^ e+) = (0.84 ± 0.21)%. To estimate 

the contribution of electrons from t decays, we consider 
the cascades B ^ t ^ e and B ^ Ds t e, with 
branching fractions taken from The rates for the 

decays B — > J/ip — > e+e" and B ■0(25') e+e" are 
also adjusted to llH. 

These irreducible background spectra are subtracted 
from the like-sign and unlike-sign spectra after correc- 
tion for electron identification efficiency. We determine 
this efficiency as a function of p* and the polar angle 9* 
using e'^e~ — > e^e~7 events and then use MC simulation 
to estimate losses in hadronic events with higher multi- 
plicities. For p* > 0.6GeV/c, the average efficiency is 
91% with an uncertainty of 1.5% estimated from the size 
of the MC correction. A summary of the yields is given 
in Table n 



TABLE I: Unlike-sign and like-sign pair yields for 0.6 < p* < 
2.5 GeV/c and their corrections with statistical and systematic 
errors. Numbers are quoted after all selection criteria. 



e e sample 



e^e* sample 



All candidates 183493 ± 434 
continuum bkgd. 22922 ± 349 
conversion, Dalitz 2978 ± 286 ± 327 
fake e 885 ± 63 ± 423 

e from same B 3200 ± 34 ± 160 



133842 ±371 
15758 ±290 
10730 ±502 ±1177 
2229 ± 182 ± 966 



e yield 

eff. corr. e yield 



153508 ± 
169654 ± 



630 ± 558 
732 ± 2235 



105126 ± 712 ± 1523 
117192 ±803 ±2510 



irreducible bkgd. 13912 ± 92 ±1341 14512 ± 97 ±2513 



corr. e yield 



155742 ± 738 ± 2606 102680 ± 809 ± 3551 



To account for B'^B^ mixing, we determine the number 
of primary electrons in the i-th p* bin from the like-sign 
and unlike-sign pairs as 



1 - foXo 
1-2/oXo 



1 - 2/oXo 



(2) 



6 



where xo = 0.186 ±0.004 is the _B"_B° mixing param- 
eter and fo = B{T{AS) -> B^B^) = 0.490 ± 0.018 
The parameter e^, is the efSciency of the additional re- 
quirement for the unUke-sign sample as defined in Eq. ^ 
The spectrum obtained from Eq.[5]is corrected for the 
effects of bremsstrahlung in the detector material using 
MC simulation. Since this correction significantly im- 
pacts the first moments, 3% for = 0.6 GeV and 0.5% 
for Eq = 1.5 GeV, we have verified that the detector ma- 
terial is simulated to better than 3%. Fig. |21 shows the 
resulting spectrum of primary electrons. 
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FIG. 2: Electron momentum spectrum from B Xev{'y) de- 
cays in the T(4S) frame after correction for efficiencies and 
bremsstrahlung, with combined statistical and systematic er- 
rors. 



Charmless semileptonic B — > X^ev decays are mod- 
eled as in |l^] by a combination of semileptonic decays 
with resonant and non-resonant hadronic systems. Us- 
ing B{B Xueu) = (2.2 ± 0.5) x lO'^ ^ to cor- 
rect for this background, we determine the moments 
Mn = YlikPk^b^cl Sfc ^b-*c whcrc k runs over all bins 
above the energy Eq and pk are the bin centers for 
n = 1 and the bin centers shifted by Mi for n — 2,3. 
These moments are then transformed into E^ moments 
M„ by correcting for the movement of the B mesons 
in the center-of-mass frame. Further biases due to the 
event selection criteria and binning are estimated from 
MC simulation. The spectra and moments presented 
are those oi B ~> XcCiyl^) decays with any number of 
photons. The moments as a function of Eq are shown 
in Fig. 13 and Table lists the principal systematic er- 
rors for Eq ~ 0.6 and 1.5 GeV. Without subtraction 
oi B ^ X^ev decays, we measure Af{'^^(1.5 GeV) — 
(1779.0 ± 1.9 ± 0.7) MeV, which is consistent with a re- 
cent measurement by CLEO Measurements with 
Eo=0 GeV have been performed by DELPHI 

We determine the partial branching fraction as 
(Efe ^bUc.u)/ {^tag eevt Ecuts), where k runs over aU bins 
withEe > Eq, Ntag = (3616.8 ± 3.5(stat.) ± 21.8(syst.)) 
xlO"^ is the background-corrected number of tag elec- 



trons, Cevt = (98.9 ±0.5)% refers to the relative efficiency 
for selecting two-electron events compared to events with 
a single etag, and Ccuts — (82.8 ± 0.3)% is the acceptance 
for the signal electron for Eq — 0.6 GeV. The result, 

B{B ^ Xev{-i), Ee > 0.6 GeV) 

= (10.36 ± 0.06(stat.) ± 0.23(syst.))%, 

is consistent with our previous measurement 01 , with the 
overall error improved by 25%. The partial branching 
fraction can be extrapolated to i?o = as part of a com- 
bined fit of the Heavy Quark EffectiveJTheory (HQET) 
parameters to the full set of moments 4] . 

Current theoretical predictions on the lepton energy 
moments do not in corp orate photon emission. Therefore 
we use PHOTOS 15] to simulate QED radiation and 
correct the moments for its impact. We verify that ra- 
diation that is not included in PHOTOS, e.g. additional 
hard photons, have no significant effect on the moments. 
The radiatively corrected moments and the estimated 
PHOTOS uncertainty ^ are given in Table |lll The 
complete listing of all moments and the full correlation 
matrix, with and without PHOTOS corrections can be 
found in Tables ITlIIVI For fitting purposes, a set of ta- 
bles and matrices with a precision of 5 significant digits 
can be obtained from the authors. 

In summary, we report a measurement of the electron 
energy spectrum of the inclusive decay B — > Xeiy and its 
branching fraction for electron energies^ above 0.6 GeV, 
which supersedes our previous result 7|. We have also 
derived branching fractions, first, second, and third mo- 
ments of electron energy spectrum from B — > XcBv de- 
cays for cut-off energies from 0.6 to 1.5 GeV. This set of 
moments combined with hadron mass moments 3] will be 
used for a significantly improved determination of HQET 
parameters and of |Vcfc| 4J. 
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FIG. 3: Measured moments of the inclusive electron energy spectrum of B — > Xccvl-y) decays as a function of the cut-off energy, 
(a) B, (b) Ml, (c) Ma and (d) M3. 



TABLE II: Results and breakdown of the systematic errors for B = tb J^^{dr/dEe) dEe , and the moments Mi, M2, and M3 
for B — > XcCV in the Z?-riioson rest frame for two A'ahu's of En. 





B[W-^] 


Mi[MeV] 


M2[10" 


-3 GeV^] 


MsflO" 


^ GeV^] 


i^olGcV] 


0.6 


1.5 


0.6 


1.5 


0.6 


1.5 


0.6 


1.5 


conversion and Dalitz pairs 


0.029 


0.001 


1.6 


0.02 


0.6 


0.00 


0.06 


0.00 


e identification efficiency 


0.151 


0.044 


2.5 


0.30 


0.6 


0.07 


0.29 


0.08 


c from same B 


0.019 


0.000 


1.3 


0.00 


0.6 


0.00 


0.03 


0.00 


B^ Ds-^e 


0.074 


0.001 


4.1 


0.04 


1.6 


0.00 


0.14 


0.00 


B^D^e 


0.060 


0.000 


3.8 


0.00 


1.6 


0.00 


0.01 


0.00 


B^T^e 


0.032 


0.002 


1.4 


0.05 


0.4 


0.00 


0.13 


0.00 


e from J/%p or '4>{2S) 


0.002 


0.001 


0.0 


0.01 


0.0 


0.01 


0.00 


0.00 


Secondary tags 


0.053 


0.011 


1.5 


0.06 


0.5 


0.00 


0.06 


0.00 


X 


0.034 


0.021 


0.8 


0.01 


0.3 


0.00 


0.03 


0.00 


tracking efficiency 


0.084 


0.033 


1.0 


0.06 


0.3 


0.02 


0.07 


0.00 


bremsstrahlung correction 


0.011 


0.028 


1.9 


0.43 


0.0 


0.05 


0.19 


0.00 


event selection 


0.052 


0.024 


0.6 


0.14 


0.0 


0.03 


0.07 


0.01 


b u subtraction 


0.047 


0.030 


1.2 


1.24 


0.6 


0.48 


0.20 


0.17 


B momentum correction 


0.000 


0.005 


0.0 


0.19 


0.1 


0.10 


0.04 


0.02 


Ntag normalization 


0.068 


0.030 














Moments 


10.17 


4.54 


1414.3 


1769.2 


148.5 


29.8 


-9.97 


2.11 


±(stat.) 


0.06 


0.03 


3.7 


1.8 


2.0 


0.8 


0.79 


0.44 


±(sys.) 


0.23 


0.08 


7.4 


1.4 


2.7 


0.5 


0.48 


0.20 


Moments with rad. correction 


10.30 


4.79 


1432.8 


1774.3 


148.0 


30.3 


-12.05 


2.12 


±(stat.) 


0.06 


0.03 


3.9 


1.9 


2.2 


0.9 


0.88 


0.47 


±(sys.) 


0.24 


0.09 


7.8 


1.4 


3.1 


0.5 


0.46 


0.20 
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TABLE III: Measured moments Mi, M2, M3, and B for five cut-off energies Eq with first their statistical error and second 
their systematic error. The moments are given in the i?-meson rest frame and are defined to include Xc hadronic states only 
and to include decays B Xcevy with any number of photons. 



So[GeV] 


B[io-^] 


Mi[MeV] 


M2[10-^GeV2] 


M3[10-^GeV^] 


0.6 


10.17 ± 0.06 ± 0.23 


1414.3 ± 3.7 ± 7.4 


148.5 ± 2.0 ± 2.7 


-9.97 ± 0.79 ± 0.48 


0.8 


9.43 ± 0.05 ± 0.19 


1469.8 ± 2.4 ± 3.6 


117.3 ± 1.2 ± 0.9 


-1.82 ± 0.60 ± 0.39 


1.0 


8.42 ± 0.04 ± 0.16 


1537.2 ± 1.9 ± 2.1 


88.3 ± 0.9 ± 0.6 


1.95 ± 0.53 ± 0.30 


1.2 


7.05 ± 0.04 ± 0.13 


1621.3 ± 1.8 ± 1.6 


61.2 ± 0.9 ± 0.6 


2.97 ± 0.49 ± 0.24 


1.5 


4.54 ± 0.03 ± 0.08 


1769.2 ± 1.8 ± 1.4 


29.8 ± 0.8 ± 0.5 


2.11 ± 0.44 ± 0.20 



TABLE IV: Measured moments Mi, M2, M-j, and B for five cut-ofi' energies Eo with first their statistical error and second their 
systematic error. The moments are given in the _B-mcson rest frame and are corrected for QED radiative effects. 



Eo[GeV] 


B[W-'^] 


Mi[MeV] 




M2[10-^GeV2] 


M3[10~^GeV3] 


0.6 


10.30 ± 0.06 ± 0.24 


1432.8 ± 3.9 ± 7.8 


148.0 ± 2.2 


± 3.1 


-12.05 ± 0.88 ± 0.46 


0.8 


9.61 ± 0.05 ± 0.20 


1484.8 ± 2.6 ± 


3.7 


117.7 ± 1.3 


± 1.0 


-3.18 ± 0.64 ± 0.38 


1.0 


8.65 ± 0.04 ± 0.17 


1548.7 ± 2.0 ± 


2.2 


89.1 ± 1.0 


± 0.6 


1.19 ± 0.57 ± 0.30 


1.2 


7.31 ± 0.04 ± 0.14 


1629.9 ± 1.9 ± 


1.7 


62.1 ± 0.9 


± 0.6 


2.66 ± 0.52 ± 0.24 


1.5 


4.79 ± 0.03 ± 0.09 


1774.3 ± 1.9 ± 


1.4 


30.3 ± 0.9 


± 0.5 


2.12 ± 0.47 ± 0.20 



TABLE V: Correlation matrix of the 20 measured moments. Matrix elements are for the full errors (statistical and systematic added in quadrature) in percent. 
The superscript i refers to the energy thresholds Ei = 0.6, 0.8, 1.0, 1.2, 1.5 GeV for i = 1, 2, 3, 4, 5 respectively. 



Ml Ml 


Ml 


Mt 


M't 


Ml 


M| 


M| 


Ml 


Ml 


Ml 


M| 


M| 


M| 


M| 




B^ 


B^ 


g4 


B« 


Ml 100.0 83.2 


64.0 


48.3 


34.0 


-73.6 


-26.5 


12.3 


22.5 


23.4 


13.5 


-9.2 


0.7 


10.5 


15.7 


-50.8 


-26.4 


-11.8 


-2.5 


8.0 


M^ 100.0 


81.4 


66.8 


51.5 


-35.1 


-20.0 


27.3 


37.8 


38.3 


-7.3 


-0.3 


9.7 


21.3 


27.2 


-36.5 


-24.1 


-7.4 


3.3 


16.7 


Mi 


100.0 


84.0 


68.5 


-3.5 


26.6 


40.4 


53.9 


54.5 


5.5 


6.6 


21.0 


34.3 


41.1 


-16.7 


-7.2 


-0.8 


12.0 


27.8 


Mt 




100.0 


79.1 


15.5 


49.0 


69.8 


61.9 


65.7 


23.8 


28.0 


33.9 


43.2 


51.6 


-0.3 


6.7 


12.1 


15.9 


35.1 


Ml 






100.0 


28.7 


63.1 


85.6 


88.8 


79.7 


48.7 


59.7 


64.4 


66.4 


65.1 


17.6 


23.4 


27.5 


31.5 


36.3 


Ml 








100.0 


66.5 


42.8 


35.2 


33.4 


3.9 


41.9 


38.0 


33.4 


30.9 


66.1 


44.1 


33.9 


29.6 


28.0 


Ml 










100.0 


75.2 


68.9 


67.3 


59.0 


61.0 


63.1 


61.4 


60.6 


53.9 


51.4 


40.3 


38.0 


41.1 


M| 












100.0 


90.9 


88.6 


66.7 


81.3 


79.7 


78.9 


78.8 


36.3 


39.4 


41.9 


39.3 


45.1 


Ml 














100.0 


94.5 


66.3 


82.1 


87.4 


87.3 


84.9 


30.3 


35.3 


38.4 


41.9 


44.0 


Ml 
















100.0 


67.4 


83.8 


90.8 


94.5 


93.6 


29.0 


34.3 


37.3 


40.1 


45.3 


Ml 


















100.0 


79.0 


76.8 


72.2 


68.2 


31.4 


43.7 


41.6 


38.2 


34.4 


Ml 




















100.0 


95.9 


90.5 


85.2 


47.5 


48.3 


48.5 


44.5 


39.5 


M| 






















100.0 


97.4 


92.8 


41.4 


44.0 


44.8 


44.0 


40.0 


M| 
























100.0 


97.5 


33.7 


37.4 


39.1 


40.6 


39.5 


Ml 


























100.0 


28.3 


32.4 


34.5 


36.1 


39.3 






























100.0 


95.4 


90.4 


86.3 


80.1 
































100.0 


97.7 


95.1 


90.1 


B^ 
































100.0 


98.3 


94.3 


B" 


































100.0 


96.5 


B' 




































100.0 
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